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ABSTRACT 

We analyze star formation (SF) as a function of stellar mass (Af*) and redshift z in the All Wave- 
length Extended Groth Strip International Survey (AEGIS). For 2905 field galaxies, complete to 
10 10 (10 10 8 )Mq at z < 0.7(1), with Keck spectroscopic redshifts out to z = 1.1, we compile SF rates 
(SFR) from emission lines, GALEX, and Spitzer MIPS 24/im photometry, optical-NIR M* measure- 
ments, and HST morphologies. Galaxies with reliable signs of SF form a distinct "main sequence 
(MS)", with a limited range of SFR at a given M* and z (la < ±0.3dex), and log(SFR) approxi- 
mately proportional to log (A/*). The range of log(SFR) remains constant to z > 1, while the MS as a 
whole moves to higher SFR as z increases. The range of SFR along the MS constrains the amplitude 
of episodic variations of SF, and the effect of mergers on SFR. Typical galaxies spend ~ 67(95)% of 
their lifetime since z — 1 within a factor of < 2(4) of their average SFR at a given M* and z. The 
dominant mode of the evolution of SF since z ~ 1 is apparently a gradual decline of the average SFR 
in most individual galaxies, not a decreasing frequency of starburst episodes, or a decreasing factor by 
which SFR are enhanced in starbursts. LIRGs at z ~ 1 seem to mostly reflect the high SFR typical 
for massive galaxies at that epoch. The smooth MS may reflect that the same set of few physical 
processes governs star formation prior to additional quenching processes. A gradual process like gas 
exhaustion may play a dominant role. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
starburst 



1. INTRODUCTION 

Deep galaxy surveys have found consistently that the 
star formation rate (SFR) per unit stellar mass (M*) 
depends strongly on both Af* and redshift z, with the 
bulk of star formation (SF) occurring earlier in massive 
galaxies than in less massive systems (e.g. Guzman et 
al. 1997, Brinchmann & Ellis 2000, Juneau et al. 2005, 
Bauer et al. 2005, Bell et al. 2005, Perez-Gonzalez et 
al. 2005, Feulner et al. 2005, Papovich et al. 2006, 
Caputi et al. 2006, Reddy et al. 2006). High-SFR ob- 
jects are observed to be more abundant at higher z; it 
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is often assumed that a part of these reflect a greater 
frequency of merger-driven starburst episodes at earlier 
times. However, a comprehensive observational picture 
of the relationship between SF and mass to z ~ 1, in- 
cluding objects with a wide range of both masses and SF 
rates, has been lacking. 

This Letter is part of a series of papers that study the 
evolution of SFR and Af* in field galaxies out to z = 1.1 
in the All Wavelength Extended Groth Strip Interna- 
tional Survey (AEGIS). We combine SFR measurements 
from deep Spitzer MIPS 24/im imaging, Keck/DEEP2 
spectra, and GALEX UV photometry, allowing us both 
to recover obscured SF in IR-luminous galaxies and to 
include lower-SFR objects not detected at 24/im. Using 
optical-NIR derived Af* measurements, we analyze the 
evolution of SFR as a function of Af* and z\ we also ana- 
lyze HST/ ACS imaging and rest-frame colors to support 
the interpretation of SFR indicators. 

We adopt a a concordance cosmology (H = 
70kms- 1 Mpc" 1 , fl M = 0.3, 51a = 0.7). Values of Af* 
and SFR are based on a Kroupa (2001) IMF, following 
recent results by Hopkins & Beacom (2006). 

2. DATA SET 

Our sample includes all field galaxies with DEEP2 
spectroscopic redshifts z < 1.1, in the area where 
Spitzer MIPS 24/tm photometry and K band imaging 
to 22ABmag are available; see Davis et al. 2006 (this 
volume). Stellar masses were obtained from SED fits to 
optical/NIR photometry by Bundy et al. (2006); errors 
are < 0.3 dcx, with a mean and rms of 0.1 and 0.05 dex 
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and < 4% of errors > 0.2 dex. Fig. 1 shows data in the 
M* range where the sample is > 80% complete, adopt- 
ing the completeness analysis by Bundy et al. (2006; see 
also Cimatti et al. 2006), for a total of 2905 galaxies. 
We draw conclusions only where the sample is > 95% 
complete (M* > 10(10.8) dex M Q for z < 0.7(1), see ver- 
tical lines in Fig. 1). For galaxies with robust 24/im 
detections (/ > 60/iJy), SFRs were derived following Le 
Floc'h et al. (2005),' using Chary & Elbaz (2001) SED 
templates; using templates from Dale & Helou (2002) 
yields no significant differences. We then add to the 
24/um-based SFR the SFR estimated from DEEP2 emis- 
sion lines (Ha, H/3, or [Oil] 3727, depending on z) with 
no extinction correction, to account for SF from unob- 
scured regions. This approach is similar to that em- 
ployed by Bell et al. 2005; utilizing rest-frame UV contin- 
uum SFRs (as they did) in place of emission-line fluxes 
yields consistent results. Galaxies below the 24/^m de- 
tection limit are not dominated by highly extincted SF; 
for these, we use extinction-corrected SFR from emission 
lines; these can probe to roughly 10 x lower SFR than 
the 24/xm data, and are slightly more sensitive at high 
z and cover a larger area than GALEX data. Emission 
line luminosities (as calculated in Weiner et al. 2006) 
were transformed to an Ha luminosity using average line 
ratios measured from DEEP2 data (H/3/Ha = 0.198, 
[OII]/Ha = 0.69; Weiner 2006, private communication), 
and transformed to SFR using the Ha calibration of Kcn- 
nicutt (1998). The DEEP2 H/3/Ha ratio corresponds to 
an extinction of 1.30 mag at Ha assuming case B recom- 
bination, which was applied to correct the emission line 
SFRs. We use fixed rather than M#-dependent line ra- 
tios (a la Weiner et al.), because these predict extinction 
corrected SFR slightly in excess of the 24/zm derived SFR 
for high mass galaxies. Our simple but robust approach 
yields results in good agreement with SFR derived from 
GALEX data, extinction-corrected based on UV spectral 
slopes. 

For objects with both < 60/iJy and emission- 

line S/N < 2, we estimate a 2-a upper limit on SFR from 
the most sensitive emission line available, by adding 2a 
to the measured uncertain SFR, or, for non-detections, 
to the limit of S/N > 2-detectable emission line SFR 
at the galaxy's redshift. We again apply Aho = 1-30 
for extinction corrections, certainly an overestimate since 
extinction is lower in more weakly SF galaxies (Hopkins 
et al. 2001). 

We have performed a suite of tests of these SF es- 
timates, finding that adopting different SFR tracers 
changes results moderately (Noeske et al 2007, in prep.); 
qualitative results are unaffected. Random errors in our 
24/im-bascd SFRs are < 0.1 dex from photometry and 
~ 0.15 dex from scatter in the f(24/im) to L(IR) conver- 
sion (see Marcillac et al. 2006), yet total random errors 
are expected to be 0.3-0.4 dex (sec Bell et al. 2005). 
For extinction-corrected emission line SFRs, random er- 
rors are ~ 0.35 dex, including scatter about the assumed 
mean extinction. 

3. RESULTS 

Fig. 1 shows SFR as a function of M* in four indepen- 
dent redshift bins. The following discussion refers only to 
the stellar mass range where the sample is > 95% com- 
plete, marked by the vertical dotted lines in each redshift 



bin. We identify three different categories of galaxies: 

(1) The majority of galaxies show clear signs of SF, ei- 
ther robust 24/zm detections, or, at lower M*, blue colors 
and emission lines (blue symbols in Fig. 1). Quantitative 
HST morphologies (Gini/M20, Lotz et al. 2006; CAS, 
Conselice 2003) classify < 25% of these galaxies as early 
types (E,S0,Sa), and > 90% show visual signs of SF such 
as blue regions and dust lanes. Most of them lie on the 
"blue cloud", (e.g. Willmer et al. 2006), though some 
of the massive ones are red, likely dusty, star-forming 
galaxies (Bell et al. 2005). This category (blue symbols 
in Fig. 1) comprises 67(56)% of the sample at z < (>)0.7 
in the M* range where the sample is complete. 

(2) Clearly separated are galaxies without robust 24/im 
(> 60/^Jy) or emission line (S/N > 2) detections (or- 
ange arrows in Fig. 1). The upper limits on their SFR 
are conservatively high (§2), such that the true sepa- 
ration between the sequence and the other galaxies is 
likely larger than it appears here. Almost all (> 95%) of 
these galaxies are on the red sequence, and > 90(80)% at 
z < (>)0.7 have early-type quantitative morphologies in- 
cluding early-type mergers, while > 90% at z > 0.7 have 
early-type visual morphologies with no hints of current 
SF. These galaxies contribute 29(30)% of the sample at 
z < (>)0.7. 

(3) Scattered below the star-forming sequence are 
galaxies with robust emission line detections but no sig- 
nificant 24/xm emission, 5(14)% of the sample at z < 
(>)0.7. All of these galaxies (green crosses in Fig. 1) 
arc on the red sequence and their Ha, H(3 emission line 
equivalent widths tend to be low (few A). Yan et al. 
(2006) and Weiner et al. (2006) showed that the bulk 
of the line emission in red galaxies out to intermediate 
redshifts is due to LINER/ AGN emission, not SF. We 
find that 75% of those galaxies with [Oil] and H/3 de- 
tections show LINER-like line ratios, and > 55(70)% at 
z < (>)0.7 have early-type quantitative and visual mor- 
phologies that are typical for local LINERs (Yan et al. 
2006). Line emission in these red galaxies thus appears to 
be dominated by LINERs/AGN, particularly at z > 0.7 
where they are more frequent. Their SFR, derived from 
emission lines (Fig. 1), will mostly be overestimated. 
However, we find visual signs of SF in the HST images 
in < 30% of these galaxies, comparable to the fraction of 
non-early quantitative morphologies; these may be dom- 
inated by SF. 

The star forming galaxies form a distinct sequence of 
SFR with M+, which we term the "main sequence" (MS). 
The red lines in Fig. 1 enclose 34% of galaxies both above 
and below the median (red points), and thus indicate 
the equivalent of ±1<7 for a Gaussian distribution. The 
width of the MS measured in this way (the range in SFR 
about the median at a given M*) is about ums = 0.35 
dex, and seems to remain approximately constant in our 
sample over the redshift range 0.20 < z < 1.1. Sub- 
tracting a lower limit of non-systematic scatter in SFR 
(~ 0.2 dex, §2) in quadrature yields an upper limit of 
~ 0.3 dex on the intrinsic scatter, which is still broadened 
by the width of the z bins, and by additional spread from 
combining different SFR tracers. Errors in M* hardly af- 
fect ctms- 

The slope of the MS is shallower than unity, 
\og(SFR) = (0.67 ± 0.08)log(M*) - (6.19 ± 0.78) for M* 
between 10 10 and lO n M , and z = 0.2 - 0.7. There is 
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Fig. 1. — SFR vs M+ for 2905 galaxies in the EGS, in the M+ range where the data are > 80% complete; see §2 . The dotted vertical 
line marks > 95% completeness. Filled blue circles: Combined SFR from MIPS 24/xm and DEEP2 emission lines. Open blue circles: No 
24/im detection, blue U — B colors, SFR from extinction-corrected emission lines. Green crosses: Same as open blue circles, but red U — B 
colors, mostly LINER/AGN candidates (§3). Orange down arrows: No robust detection of /(24/mi) or emission lines; conservative SFR 
upper limits shown. There is a distinct sequence formed by fiducial SF galaxies (open and filled circles); galaxies with little or no SF lie 
below this sequence. Red circles show the median of log(SFR) in mass bins of 0.15 dex for main sequence galaxies (blue circles). Red 
lines include 34% of the main sequence galaxies above and 34% below the median of log(SFR), ±ltr in the case of a normal distribution. 
Horizontal black dashed line: SFR corresponding to the 24/im 80% completeness limit at the center of each z bin. 24/tm-detected galaxies 
above the magenta dot-dashed line are LIRGs (§4.2). 



a trend for the slope to flatten to higher z, but the com- 
pleteness limits do not allow a robust quantification. A 
further important result is that the normalization of the 
main sequence evolves strongly over the redshift range of 
our sample; the median SFR at fixed M* evolves down- 
wards by a factor of 3, measured at IO^Mq, from the 
our highest (median z — 0.98) to our lowest (median 
z = 0.36) redshift bin. Importantly, it appears that the 
whole of the main sequence shifts downwards with time, 
rather than just the upper envelope decreasing, which 
was also reported by a recent GALEX study at z = 0.7 
(Zamojski et al. 2007). A straightforward interpretation 
of these observations is that normal star-forming galax- 
ies possess a limited range of SFR at a given M± and z, 
which is presumably set by whatever physical processes 
regulate SF in quiescent disks. Galaxies that are not on 
the main sequence, in categories (2) and (3) above, are 
observed during or after quenching of the SF activity, 
with either low-level or no current SF, or LINER/AGN 
activity. 

4. DISCUSSION 

4.1. Completeness: Is the Main Sequence Real? 

It is obviously crucial to determine whether the "main 
sequence" that we have identified is real, or could be 
caused by selection effects or observational biases. We 
address the following possible causes of incompleteness 
or bias in our sample, again restricting the discussion to 
the range where we claim that the sample is > 95% 
complete: 

(1) Could the optically selected DEEP2 parent sample 
be missing a significant number of galaxies, or are there 
galaxies in the DEEP2 sample that lack a successful red- 
shift determination because of low S/N? (2) Could we 
be significantly underestimating the SFR in galaxies in 
our sample due to biases in our SF indicators? 

(1) The DEEP2 spectroscopic selection (R AB < 24.1) 
has been shown to be complete in the M* ranges in- 
dicated in Fig. 1 (vertical lines), from comparisons to 
various surveys with spectroscopic and deep photometric 
redshifts, including in particular the -fT-selected K20 sur- 



vey, which should be less affected by extinction (Willmer 
et al. 2006, Bundy et al. 2006, Cimatti et al. 2006). 
For galaxies that are below our 24/mi detection limit, we 
expect the extinction to be moderate, and would expect 
these galaxies to be picked up in if-selected surveys, but 
no such population is found to be missed by DEEP2. 

More obscured populations can be probed through the 
deep Spitzer IRAC 3.6/im data in AEGIS, which at a 
given redshift are a proxy for M+ yet are barely affected 
by extinction. We have compared the distribution of 
/(24/xm) at a given /(3.6^tm) and z in the DEEP2 Rab~ 
selected sample and an IRAC /(3.6/^m)-selected sample 
with IRAC-based photometric redshifts. We find no ev- 
idence that DEEP2 misses a significant population of 
heavily obscured, star- forming galaxies at z < 1, which 
could populate the area above the upper boundary of the 
MS. This agrees with the results of Houck et al. (2005) 
and Weedman et al. (2006) in the large area NDWFS, 
which indicate that such missed /(24^m)-bright, opti- 
cally faint galaxies at z < 1 would contribute < 1% of 
our sample. 

(2) The 24/im completeness limit (horizontal black 
dashed line in Fig. 1) intersects the main sequence 
in each redshift bin. As discussed in Sj3l most galax- 
ies below the MS are red, early type, non-SF, and/or 
LINER/AGN dominated (shown as orange arrows and 
green crosses). However, a fraction show spiral/late- type 
morphologies or visual signs of possible SF (§3). In prin- 
ciple, these red galaxies could have dust-obscured SF, 
unrecovered by emission lines, yet lie below the 24/^m de- 
tection limit. Their true SFRs could then be anywhere 
up to the 24/im limit, in which case they may not be a 
distinct population, but rather a downward continuation 
of the MS. If this were the case, these galaxies would 
make up < 10(20)% of the MS at z < (>)0.7. We can 
constrain the maximal effect of missed, dust-obscured SF 
in these galaxies on the ler range of SFR along the MS 
by including in the calculation of o~ms all red, 24/xm- 
undetected galaxies with spiral/late- type morphologies, 
and Ha,H/3, and/or [Oil] line emission down to spurious 
detections (i.e., 100% error in EW). For the extremes of 
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either only the emission line SFR or the maximal SFR, 
corresponding to the 24/im limit, the measured width of 
the MS increases by ~ 0.05 dex or not at all, respectively. 

Thus we argue that the relatively sharp upper limit of 
the MS is real, as our selection does not miss obscured 
sources with high SFR. The sharpness of the lower limit 
is more uncertain with our current data, but we find 
that only a small fraction of galaxies that we placed be- 
low the MS could have underestimated SF rates which 
would "blur out" the lower edge of the MS. Very deep 
24^ data at z ~ 1 from GOODS (Elbaz 2006, private 
communication) unambiguously confirm this result, par- 
ticularly a well-defined lower boundary to the sequence. 

4.2. Constraints on Episodic Star Formation 

Studies based on local samples (Brinchmann et al. 
2004, Salim et al. 2005 for SDSS; Lee 2006) have il- 
lustrated a relationship between SFR and M*, and have 
identified two populations: galaxies on a star-forming 
sequence, and "quenched" galaxies with little or no de- 
tectable SF. At higher z, previous studies (see §1) had 
merely described an upper envelope of SFR in the SFR- 
M* diagrams. We have employed a variety of SFR trac- 
ers and other evidence from AEGIS to show that the 
SF sequence persists out to z ~ 1, with a similar dis- 
persion in log(SFR) at fixed M* but with a decrease in 
normalization of a factor of 3 from z = 0.98 to z = 0.36, 
measured at M* = 10 11 M Q . The global star formation 
rate density has also decreased by a factor of 3 over this 
same interval (Hopkins et al. 2004, Eq. 3). One possi- 
ble physical explanation for this decline is a decreasing 
contribution from starbursts in gas-rich galaxy mergers. 
However, if this were the dominant factor causing the 
decline, we would expect to see the upper envelope of 
the main sequence move downwards with time, with the 
region populated by "normal" galaxies maintaining the 
same normalization. This is contrary to what we see 
in AEGIS: the region of the SFR-M* space populated 
by main sequence galaxies at z ~ (Brinchmann et al. 
2004) is nearly empty at z ~ 0.7-1, though these galaxies 
should be detectable in our survey. 

We can use our observed MS to quantitatively con- 
strain the duty cycle of episodic variations of SFR around 
an average level. We adopt the densely populated peak 
of the SFR distribution (the median) as this baseline 
level. The 1(2)<7 ranges about the median of log(SFR), 
±0.3(0.6) dex, include 68(95)% of the galaxies. We can 
hence infer that SFR variations exceeding ±0.3(0.6) dex, 
factors of 2(4), have duty cycles < 32(5)%. These cor- 
respond to total times of < 2.5(0.4) Gyr since z = 1. 
The amplitude of these variations (< x4) is consistent 
with gas-poor or minor mergers, rather than the peak 
SFR of gas-rich major mergers (Springel 2000, Cox et 
al. 2006). Excursions in SFR > 5x above the median 
are rare, ~ 1%, consistent with galaxies spending < 100 
Myr in such strong burst episodes since z = 1. Of course 
these arguments are only valid in a statistical sense: a 



fraction of galaxies could have a lower average level of 
SF and undergo larger excursions, but at the expense of 
reducing the allowed range of SFR for the remainder of 
the population. 

Previous studies have found that < 30% of SF at 
z ~ 0.7 occurs in morphologically disturbed galaxies 
(Wolf et al. 2005; Bell et al. 2005) or close pairs (Lin et 
al. 2006). Semi-analytic models predict that about 5% of 
the SF at z ~ 0.7 is due to major mergers, with the con- 
tribution due to minor mergers being more uncertain, 
but ranging from ~ 11-45 % (Somcrville et al. 2001; 
Wolf et al. 2005). These direct constraints and theoreti- 
cal expectations are consistent with the conclusions that 
we have drawn here from the SFR-M* distributions. A 
related comment pertains to the nature of Luminous In- 
frared Galaxies (LIRGs, £(8-1000/x) > 10 n i Q ). LIRGs 
at z ~ are rare, mostly interacting galaxies (Sanders & 
Mirabel 1996) with strong starbursts (SFR > 5x above 
those of typical spirals). At z ~ 1, LIRGs seem to mostly 
represent the high level of SFR in almost all massive SF 
galaxies, rather then extreme starbursts (Fig. 1) . 

In summary, we suggest a picture in which we are wit- 
nessing a gradual decline in the SFR of most galaxies 
since z ~ 1, accompanied by rapid quenching in a frac- 
tion of (massive) galaxies. Presumably the regularity and 
constant dispersion of the main sequence out to z ~ 1 
means that the same physics that regulates SF in local 
disk galaxies is operating, indicating significant evolu- 
tion either in the gas supply or SF efficiencies over this 
interval. 

In the accompanying Letter (Noeske et al. 2007, this 
volume) , we show that the slope and evolution of the MS 
can be understood as gradual gas exhaustion in a model 
in which galaxy age and SF timescale are a function of 
galaxy mass, and the dispersion of the MS is interpreted 
as resulting from a spread in age and SF timescales at a 
given mass. 
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